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First Simulation Results for the Photon in a Non-Commutative Space 
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We present preliminary simulation results for QED in a non-commutative 4d space-time, which is discretized 
to a fuzzy lattice. Its numerical treatment becomes feasible after its mapping onto a dimensionally reduced 
twisted Eguchi-Kawai matrix model. In this formulation we investigate the Wilson loops and in particular the 
Creutz ratios. This is an ongoing project which aims at non-perturbative predictions for the photon, which can 
be confronted with phenomenology in order to verify the possible existence of non-commutativity in nature. 
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1. THE NON-COMMUTATIVE U(l) 
MODEL 

We investigate non-perturbatively QED on a 
non-commutative (NC) space-time where a 
commutation relation of the form 



(1) 



is introduced. This has the effect of replacing 
points in regular space with "Planck cells" . O^j/ is 
a real and antisymmetric non-commutativity ten- 
sor. 

Here we address pure U(l) gauge theory in a 
Euclidean NC space-time. Our main tool is the 
star product, which is defined as 

f{x) ★ g{x) = f{x) exp (^^% ©p'- 9{x) ■ (2) 

This definition absorbs the "space deformation" 
into the products of fields, thereby allowing us 
to treat the space-time operators x and functions 
thereof as regular space-time coordinates again, 

4){x)ip{x) • (3) 

Equipped with this tool we are able to write 
down the gauge action of the pure U(l) theory as 

^[^] = -i^/ d*xTT[F^,{x)*F^,ix)] , (4) 
where F^^, is given by 

F^,^ = d^A^ - d^A^ - i{A^, ^A^-A^-k A^). (5) 



The gauge fields in this theory do not commute, 
as the Yang-Mills term shows. Another unusual 
feature of this theory are complex, ★-gauge in- 
variant Wilson loops [2] • Moreover a 0-deformed 
dispersion relation for the photon is expected [H]. 
Such a modification was observed for the NC Ac/)** 
model Uj . This is due to a remarkable feature of 
the NC field theories, the UV/IR mixing of diver- 
gences IH] ■ That property makes the perturbative 
treatment beyond one loop extremely difficult. 

2. LATTICE FORMULATION OF THE 
NC PLANE 

We carried out simulations involving two NC 
directions, as well as two commutative directions, 
which include the Euclidean time. To be explicit, 
we set Oi2 = —621 — 9, and ~ otherwise. 
The parameter 9 was assumed to be constant 
throughout the space-time. 

In order to perform Monte Carlo simulations 
one needs to formulate the theory on a lattice. 
This is achieved by the operator relation |H] 



exp [ a;^ ) = 1, 



= 1,2 , 



(6) 



which imposes a discrete non-commutative space- 
time. After having obtained the lattice space- 
time, it needs to be populated by U(l) link vari- 
ables. The key to render their simulation — in- 
volving the non-local ^-product — feasible is the 
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Morita equivalence of NC U(7V) gauge theo- 
ries to the twisted Eguchi-Kawai (TEK) model 
IH], with the action 



(7) 



Here the Ufj, are U(A^) matrices. The twist 



^uf, is given by exp 



N 



with n^i, G Z. In 
-ni2 = (N + l)/2 



our mapping we chose n2i 
for some odd N, and n^^ = otherwise. 

We discretized the commutative plane to a reg- 
ular (TV— 1) X (A^— 1) lattice, where on each lattice 
site the NC plane is represented by a x iV TEK 
model (A'' is odd). 

The non-commutativity parameter 9 and the 
lattice parameters TV and a are related as [H] 



n 



(8) 



In the N ^ oo, a ^ limit this theory describes 
a continuum gauge theory; which theory one ob- 
tains depends on the way this limit is taken. 
Keeping 6 finite in the continuum and infinite vol- 
ume corresponds to the double scaling limit. In 
the case of finite N and a the theory maps exactly 
to a NC U(l) lattice theory 0. 

In order to be able to apply the heat bath al- 
gorithm, the action was linearized in [/^, as a 
generalization of the scheme of Ref . [5] ■ 

3. RESULTS 

In the simulation results presented here we set 
N = 25. The upper part of Fig. ^ shows the real 
part of the plaquette values at different /3. In this 
case all six possible orientations of the plaquettes 
are averaged. The plot shows that the results 
match the anticipated asymptotic behavior 
The phase transition occurs at around P ~ 0.35. 

The decay of Wilson loops observed from dif- 
ferently sized squares is shown in the lower part 
of Fig. n In that plot the results are split ac- 
cording to the plaquette orientations (NC, mixed 
and commutative) . The real part of the NC Wil- 
son loops oscillates, while eventually decaying to- 
wards zero. This is in qualitative agreement with 
the observations in Ref. JHI for d — 2. For the 
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Figure 1. Above: The real part of the plaque- 
ttes, averaged over all possible orientations. The 
dashed lines show the asymptotic behavior ex- 
pected based on a strong/weak coupling expansion. 
Below: The real part of the plaquettes, split ac- 
cording to their orientation. 



other two types we see a monotonous decay to- 
wards zero. Notice that these results were ob- 
tained in the weak coupling phase. 

We have calculated Creutz ratios x for various 
rectangular Wilson loops. The results are shown 
in Fig. 121 From these we see that the string ten- 
sion seems to approach zero as /3 — > oo, which 
corresponds to the continuum limit. 

Fig-Elshows the correlation function of the pla- 
quettes lying completely in the NC plane, and 
separated by At in Euclidean time. These data 
were taken at /3 = 0.4, so we are barely in the 
weak coupling phase. The decay seems to be ex- 
ponential, but more statistics are required to con- 
firm this behavior. 

Fig. 01 shows the commutative analog to the 
correlation function discussed above. However, 
in that case the decay does not seem exponential. 
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Figure 2. Different size Creutz ratios x of rectan- 
gular Wilson loops, averaged over all planes, at 
various couplings. 
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Figure 3. Temporal correlation of the plaquettes 
lying completely in the NC plane. 

4. SUMMARY AND CONCLUSION 

We have simulated a NC U(l) field theory. 
First observables such as the action, plaquettes, 
Wilson loops and Creutz ratios could be mea- 
sured. The results for the action and plaquettes 
agree with asymptotic predictions. There seems 
to be a vanishing string tension at weak coupling 
in the part of the parameter space we explored. 

Further work will hopefully allow us to identify 
a physical scale and to extrapolate to a NC con- 
tinuum limit by means of double scaling (i.e. a 
fixed product Na^). Thus we hope to ghmpse at 
the dispersion relations for the NC photon, which 
might allow us to confront the theory with nature. 
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Figure 4. Temporal correlation of the plaquettes 
lying completely in the commutative plane. 

the "Norddeutscher Verbund fiir Hoch- und 
Hochstleistungsrechnen" (HLRN). 
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